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I.  Introduction 

The  long  term  objectives  of  this  research  are  to  provide  a  tool  for  observing 
oceanic  events  on  a  length  scale  smaller  than  those  associated  with  internal  waves 
and  to  provide  improved  predictive  capabilities  for  direct  path  propagation  fluctuations 
in  the  20  to  70  KHz  regime,  where  small  scale  events  can  dominate  forward  acoustic 
scattering. 

Scattering  of  an  acoustic  signal  by  the  dynamic  processes  and  the 
heterogeneous  chemical  and  physical  distribution  of  properties  that  exist  in  the  open 
ocean  have  received  considerable  attention  in  the  last  decade.  Our  understanding  of 
the  effects  of  stochastic  volume  scatter  mechanisms  on  a  propagating  acoustic  field, 
both  at  the  level  of  the  second-order  statistics  (signal  coherence  measures,  spread 
functions)  and  at  the  level  of  fourth-order  statistics  (signal  intensity  fluctuations  and 
correlations)  is  reasonably  complete.  A  comparably  detailed  description  of  the 
structure  and  mechanisms  of  physical  oceanography  at  intermediate-to-fine  scales  is 
not  yet  well  developed. 

This  report  represents  a  part  of  an  effort  to  develop  theory  on  which  to  base  the 
development  of  new  techniques  to  be  used  in  developing  a  better  understanding  of 
the  intermediate  and  fine  scale  physical  processes.  These  are  the  processes  that 
exist  between  the  spatial  scales  of  internal  waves  and  the  dissipative  range  of 
turbulent  ocean  structure.  /The  development  of  a  new  technique  or  tool  to  study  such 
processes  and  distributions  of  ocean  properties  must  include  the  validation  of  the 
mathematical  theory  on  which  the  technique  is  to  be  based.  Consequently,  some 
initial  efforts  directed  toward  the  development  of  the  procedures  which  are  necessary 
to  accomplish  that  validation  are  also  described  in  this  report. 

The  obective  of  this  research  is  to  provide  a  tool  that  has  applicability  in  many 
ocean  environments.  As  one  attempts  to  develop  a  new  tool,  it  is  often  advantageous 
to  seek  out  an  environment  that  initially  maximizes  one’s  signal  and  minimizes  the 
measurement  noise.  That  environment  should,  however,  be  representative  of  the 
ocean.  Internal  waves  often  dominate  the  energy  spectrum  of  processes  that  influence 
fluctuations  in  sound  propagation  in  the  temperate  oceans.  This  can  make  it  difficult  to 
separate  internal  wave  induced  fluctuations  from  those  caused  by  simple  stochastic 
heterogeneity  in  sound  speed.  The  motions  of  even  a  relatively  large  research  ship  in 
normal  seas  would  increase  the  difficulties  one  might  expect  to  encounter  in  the 
accurate  measurement  of  both  the  acoustic  fluctuations  and  the  oceanographic 
structure. 

We  have  suggested  the  Arctic  as  a  possible  location  in  which  such  a  technique 
as  is  proposed  could  be  tested  and  developed.  The  Arctic  has  the  advantage  that  the 
internal  wave  energy  appears  to  be  lower  than  in  ice  .ree  waters.  The  technical 
advantage  of  making  the  measurements  from  a  stable  platform  (the  ice)  is  considered 
to  outweigh  the  logistical  difficulties  of  working  in  the  Arctic  environment.  This 
suggestion  is  also  influenced  by  a  lack  of  knowledge  in  the  Arctic  about  the 
parameters  that  we  would  be  measuring.  Thus,  if  the  Arctic  were  used  as  a  natural 
laboratory  for  a  part  of  the  proposed  development,  there  would  be  an  added 
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advantage  of  potentially  gaining  direct,  new  knowledge  about  the  physical  and 
acoustical  oceanography  of  this  unique  environment. 

In  simple  terms,  our  approach  is  to  develop  those  stochastic  models  and 
computer  codes  that  describe  the  forward  scattering  process.  Those  models  will 
involve  the  parameters  of  physical  oceanography  that  we  ultimately  wish  to  measure. 
The  validation  of  the  models  requires  measurement  of  acoustic  intensity  after  it 
propagates  through  a  medium  with  a  sound  speed  heterogeneity.  It  also  requires  a 
direct  measurement  of  the  stochastic  descriptors  of  the  sound  speed  field.  After  the 
validation  of  the  "forward"  model,  the  next  step  is  the  solution  of  the  "inverse"  problem, 

i.e.,  expression  of  the  values  for  the  parameters  of  physical  oceanography  in  terms  of 
the  acoustic  measurements.  Finally,  this  "inverse"  solution  must  be  successfully 
compared  to  experimental  measurements. 

The  frequency  or  frequencies  of  the  sound  used  in  any  tool  developed  to  study 
the  physical  oceanography  of  the  ocean  must  be  carefully  selected  to  be  optimally 
sensitive  to  those  phenomena  under  investigation.  As  described  later  in  this  report,  for 
our  application,  those  frequencies  are  in  the  hundreds  of  kHz  to  the  low  MHz  range. 
Ultimately,  however,  if  the  proposed  tool  is  successful  in  developing  an  accurate 
description  of  the  basic  distribution  of  ocean  properties,  one  would  wish  to  use  the 
results  for  predictive  purposes.  Consequently,  the  final  test  we  would  envision,  would 
be  the  measurement  of  forward  scattering  at  frequencies  in  the  20-70  KHz  range  and 
comparison  of  those  measured  results  with  the  propagation  fluctuations  predicted  from 
an  estimation  of  the  ocean  properties  made  with  the  proposed  tool  (operating  at  the 
higher  frequencies).  This  would  complete  the  long  range  program  that  we  have 
envisioned. 

As  an  initial  effort  toward  the  attainment  of  the  long  term  goals  outlined  above,  we 
have  conducted  a  number  of  theoretical  and  experimental  studies  during  the  term  of 
this  contract.  These  studies  have  included  the  following  topics.  Several  additional 
topics  have  been  covered  in  an  interim  report. 

1 .  Extension  of  the  isotropic  theory  to  include  anisotropic  effects. 

2.  An  assessment  of  the  Arctic  as  a  potential  natural  laboratory  for  the  testing  of 
the  theory. 

3.  An  assessment  of  the  effect  of  random  velocity  fluctuations  on  underwater 
acoustic  scattering. 

4.  Determination  of  the  intensity  distribution  from  a  circular  source  in  an 
isotropic  random  medium. 

5.  A  study  of  the  optimum  parameters  for  a  laboratory  tank  experiment. 

6.  A  brief  description  of  the  tank  experiment. 

7.  An  analysis  of  results  of  the  tank  experiment. 
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II.  Extension  of  the  Isotropic  Theory  to  Include  Anisotropic  Effects 

During  the  course  of  this  project,  the  inverse  theory  for  isotropic  scattering  was 
extended  to  include  anisotropic  effects.  We  first  derived  the  relevant  equations  and 
direct  solutions  and  then  found  expressions  for  the  inverse  relations.  Computer  codes 
were  written  for  the  solutions  in  the  case  of  an  exponential  correlation  function.  In  the 
case  of  strong  scattering,  a  simple  analytic  solution  for  the  Gaussian  beam  was 
developed.  In  addition  to  considering  the  measurement  of  the  intensity  distribution 
associated  with  transmission  of  a  Gaussian  sound  beam,  we  also  considered 
measurement  of  the  coherence  functions  which  result  from  a  point  source.  A  computer 
code  was  written,  and  an  analytic  solution  developed,  for  the  case  of  an  exponential 
correlation  function  for  strong  scattering.  We  also  provide  analytic  expressions  for  the 
case  of  other  wave-number  power  spectra  when  the  scattering  is  strong. 

A.  Definitions 

To  treat  the  problem  from  a  stochastic  point  of  view  we  define  the  following 
coherence  functions: 


{r(xT1lxT2,z)}  =  {u(xT1,z)u*(xT2,z)}, 


{U(XT1,XT2,XT3,XT4,Z)}  = 

{u(xT1,z)u*(xT2,z)u*(xT3,z)u{xT4,z)}. 


Here,  u*  denotes  the  complex  conjugate  of  the  acoustic  field  u,  and  the  transverse 
vector  xT  has  the  components  x  and  y.  The  brackets,  {  },  denote  an  ensemble 
average. 

The  intensity,  {l(xT,z)},  is  found  from  {r(xT1  ,xT2,z)}  by  setting  xT1  =  xT2.  The  intensity 
fluctuations  {l2(xT,z)}  and  the  coherence  of  intensities  {R|(xT1,xT2,z)}  are  found  from  {L1} 
from  using  the  following  relations: 


{|2(xt,z)}  =  {U(xT,xT,xT,xT,z)}, 


{R|(Xti.XT2,Z)}  =  {LMXtlXjlX^.Xts.Z)}. 


We  write  p(x)  in  the  form  p(x)  =  p(x)  +  p'(x)  where  p(x)  represents  the  mean 
background  profile  and  p'(x)  represents  the  random  index  of  refraction  fluctuations. 
We  shall  assume  here  that  the  statistics  of  the  index-of-refraction  field  are 
homogeneous.  The  two-point  correlation  function,  a(x1tx2),  is  defined  as, 


a(x1,x2)  =  {p’(x1)p’(x2)} 
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and  depends  upon  r  =  -  x2.  In  addition,  we  define  the  function 

oo 

®(XT1  '  XT2)  =  VU)  I  °(XT1  '  xT2'rz)drz- 

-OO 


B.  Basic  Equations 

The  moment  equations  for  {T}  and  {L1}  in  the  anisotropic  case  are 
3{r}/az  =  (i/2k)[V2Tl  -  72T2Kn  - 
k2[a(0)  -  g(xT1  -  xT2)]{r}  + 

ik[fl(xT1)-fl(xT2)Kn,  (1) 

a{Ln/az  =  (i/2k)[v2Ti  -vhz  +V2T4]{L1}  + 
k2[  0(XT1  -  XT2)  +  o(XT1  -  Xj3)  + 

a(xT2  -  XT4)  +  0(X-f3  -  xT4)  - 

<j(xT1  *  XT4)  ■  s(xT2  ’  XT3)  *  (2) 

2  5  (0)]{U}  + 

•k[  p(xT1)-p(xT2)-p(xT3)  + 

ft(xT4)]{l-1}- 


These  equations  are  a  direct  generalization  of  the  isotropic  case  in  which  the  scalar 
separation  distance  r  is  replaced  by  the  vector  separation  distance, 
rij  -  xTi  *  XT|- 

We  have  included  in  the  above  equations  a  mean  background  effect  composed 
of  a  sum  of  ji(xT1)  terms.  It  can  easily  be  shown  that,  for  a  linear  profile,  this  term  has 
no  real  effect  on  our  results.  The  measured  patterns  are  simply  displaced  in  position 
or  angle. 

For  both  the  second-  and  fourth-order  coherence  functions,  it  is  convenient  to  use  a 
different  set  of  independent  coordinates.  For  the  second-order  case  we  transform  to 
sum  and  difference  coordinates  defined  as 
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p  =  (Xfi  +  xT2)/2. 

s  =  (xT1  -  XT2). 

In  terms  of  s  and  p,  Eq-  0)  becomes 
d{T)ldz  =  (i/k)  3/3 S  *  did p{H  '  k2[5(0)  -  5(S )]{0 

+  ik[p(p  +  s/2)  -  fl(p  *  s/2)]{r>- 

For  the  fourth-order  case,  we  transform  to  the  coordinates 
s  =  (1'2)(Xti  -  XT2  +  x73  -  XT4), 

P  =  (V2)(Xti  +X-T2-XT3-XT4), 
q  =  (Xti  '  Xt2  •  XT3  +  XT4)' 

R  =  (V4)(XT1  +  XT2  +  XT3  +  XT4)- 

Here, 

r12  =  s  +  q/2  ri3  -  P  +  Q12 

r14  =  s  +  p  r23  =  -s  +  p 

r24  =  p-q/2  =  s  -  q/2. 


Using  Eqs.  (5)  and  (6),  Eq.  (2)  may  be  written  as 

3{U)/3z  =  (i/k)  [(3/3s)  •  (did p)  + 

(3/3R)  ♦  (3/3q)]{L)}  +k2[  o(r12)  +  5(r13)  +o(r24)  +  rt(r34) 

a(r i4)  -  o(r23)  -  2o(0)]{L1}  + 

iK(p(xTi)  -  P-(XT2)  ‘  WXT3)  +  P-(xT4)Kl1}- 
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C.  Solutions 

1 .  Plane  Wave  Boundary  Condition:  When  there  is  no  mean  profile  (i.e.,  (I  =  0),  Eq.  (4) 
may  be  readily  solved  for  an  initial  plane  wave.  In  the  plane  wave  case  ,  the  term 

(d/dp)  {r}  =0  and  the  solution  is 

{r(s,z)}  =  I  exp  (-k2[o(0)  -  5(s)]z), 

where  I  is  the  initial  wave  intensity.  [If  there  is  a  linear  profile,  the  right-hand  side  of 
Eq.  (8)  is  multiplied  by  a  phase  term  of  the  form  exp  (ikAsz),  where  the  constant  A  is 
determined  by  the  slope  of  the  profile.] 

The  intensity  is  independent  of  z  and  this  is  readily  confirmed  by  setting  s  =  0.  As 
Isl  approaches  infinity,  we  find 

{r(oo,z)}  =  I  exp  (-k2  o(0)z).  (9) 

We  may  then  interpret  {F(<»,z)}  as  the  intensity  of  the  coherent  field  that  remains  after 
scattering.  The  quantity  Lf  =  1/[k2o(0)j  is  the  characteristic  distance  over  which 
signficant  scattering  occurs. 

The  perturbation  solution  for  {R,(xT1.xT2,z)}  obtained  from  Eq.  (7)  is 


oo 

{R|(xT1  iXT2,z)}  =  l2[l  +  (k2z/2)  If  §(k12)  exp(ik12  •  r12) 

-oo 


(1  -  (k/k122z)  sin[k122z/k])  dk12]  (10) 


where 


OO 

^(^12)  =  (^2n)2  11  °(r12)  exP(iki2  *  r12)drl2. 

—  OO 


A  linear  mean  profile  does  not  change  the  nature  of  this  solution. 

2.  Point  Source  Boundary  Condition:  The  solution  for  {r}  a  point  source  is 

1 

(r(s,z)>  =  (lz/z2)  exp(-k2z[c(0)  -  j  o(as)da]). 

0 


(11) 
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Here  lz  is  a  normalization  constant.  Further  discussion  of  this  boundary  condition  will 
be  given  following  the  determination  of  the  intensity  distribution  for  a  Gaussian  beam 
boundary  condition. 


3.  Gaussian  Beam  Boundary  Condition:  For  jl  =  0,  Eq.  (4)  may  be  solved  analytically 
for  a  Gaussian  beam  boundary  condition  For  this  case  we  have 


{r(xT1,xT2,0)}  =  I  exp[-(|xT1|2  +  |XT2|2)/2b2] 
=  I  exp[-(|p|2  +  |s|2/4)/b2] 
The  intensity,  {l(p,z)},  is  found  to  be 


{l(p,z)}  =  (I  b2/4jr)  \\  dw  exp(-iw  •  p) 

-oo 

exp[-w2b2(1  +  z2/k2b4)/4  +  Q(w,z)] 

where 

wz/k 

Q(w,z)  =  (k3/w)  |  [  a(wxt/w,wyt/w)  -  a(0,0)]dt. 

0 


(12) 


(13) 


(14) 


The  only  effect  that  a  linear  mean  profile  will  have  is  to  shift  the  maximum  value  of  the 
intensity. 


D.  Computer  Codes  for  the  Direct  Problem 

We  have  written  computer  codes  to  evaluate  Eq.  (13)  for  the  general  anisotropic 
case  and  in  the  isotropic  case.  We  have  only  to  specify  the  function  a  (r,j)  and  the 
parameters  k  and  b  in  order  to  find  the  intensity  (l(p,z)} ,  as  a  function  of  separation 
distance  p,  at  range  z.  In  Figs.  1  and  2  we  show  graphs  of  intensity  vs.  separation 
distance  for  k  =  420  nr1  (100  Hz)  and  b  =  0.34  m.  The  range  is  500  m.  The  scattering 
function  a  is  chosen  to  be  of  the  exponential  form 


a(wx,wy)  =  (p’2}  L,<  exp[-((l  sxl  /Lx)  +  (I  syl  /Ly))J. 


(15) 


where  we  have  chosen  (p'2}  =  1 0'9.  the  characteristic  horizontal  correlation  length, 
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Lx  =  100  m  and  the  characteristic  vertical  correlation  length,  Ly  =  1  m.  The  above 
numbers  appear  to  be  typical  values  for  the  fine  structure,  but  are  just  meant  here  as 
an  example.  The  exponential  form  was  chosen  somewhat  arbitrarily  since  we  could 
obtain  simplified  analytic  solutions  (presented  below)  to  check  the  numerical  program. 
We  also  performed  some  calculations  using  a  Gaussian  correlation  function  to  further 
check  the  numerical  code. 

In  Fig.  1  we  have  chosen  px  =  0  (to  obtain  a  vertical  intensity  profile)  and  in  Fig.  2 
we  took  py  =  0  (to  obtain  a  horizontal  intensity  profile).  As  a  result  of  the  high 
anisotropy,  the  two  profiles  are  very  different.  In  fact,  for  the  numbers  chosen,  the 
vertical  profile  is  principally  the  result  of  scattering  and  the  horizontal  profile  results 
principally  from  diffraction  of  the  source.  Note  that  the  abcissas  in  Figs.  1  and  2  are 
different. 


Fig.  1 :  Vertical  Intensity  Profile  - 1  (py,z)  vs  py  (px=0,  k=420  nr1, 
MOO  kHz,  b=0.34m,  {p12}=10'9,  z=500  m) 
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Px(m) 


Fig.  2:  Horizontal  Intensity  Profile  - 1  (px,z)  vs  px  (py=0,  k=420  nr1, 
f=1 00  kHz,  b=0.34m,  {p12}=10-9,  z=500  m) 


E.  Analytic  Approximation 

The  mean  free  path  for  scattering,  Lt  is  given  by  the  expression 
L,  =  1/({n*2}k2Lx 

When  the  condition  z/L,  »  1  is  satisfied,  then  the  integral  expression  for  Q(w,z),  given 
in  Eq.  (14),  can  often  be  simplified.  For  the  exponential  correlation  function  (Eq.  (15)), 
we  find  after  some  manipulation 

Q(w,z)  =  (1/2){p'2}kz2[l  wx|  +  Lxl  Wyl  /Ly],  (16) 

Using  Eq.  (16),  the  expression  for  (l(R,z)}  can  be  evaluated  in  terms  of  error 
functions.  However,  because  of  the  condition  L^Ly  »  1 ,  we  find  that, for  an  under-ice 
Arctic  experiment,  we  would  anticipate  diffraction  of  the  source  to  dominate  the 
intensity  behavior  in  the  horizontal  (x)  direction  and  scattering  to  dominate  in  the 
vertical  (y)  direction.  Under  these  conditions  we  find,  finally, 

(l(R,z)}/l  =  [kb3Nn  z]  [Kz2/(K2z4  +  Ry2)]  exp[-(Rxkb/z)2] ,  (1 7) 

where  K  =  (V2){p’2}kLx/Ly  and  we  have  assumed  that  z/kb2  »  1 . 

In  order  to  quickly  assess  the  effect  of  changing  frequency  and  range  we  can  use  Eq. 
(17).  It  must  be  emphasized,  however,  that  the  function  (l(R,z)}  is  very  sensitive  to 
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the  form  of  5.  Both  the  intensity  at  R  =  0  and  the  characteristic  decay  length  of 
theintensity  in  the  vertical  direction  may  change  appreciably  when  the  form  of  a  is 
changed. 

For  the  values  {p'2}  =  4{c'2}/{c}2  =  1 0*9,  Lx  =  100  m  and  Ly  =  1  m  we  present  a 
small  table  for  the  characteristic  vertical  decay,  Ryc,  when  Rx  =  0.  Here  z/Lf  =  8.82.  The 
value  of  Ryc  is  taken  at  the  point  at  which  the  intensity  is  half  its  initial  value.  From  Eq. 
(17),  we  find, 


Ryc  — 

{p*2}  kLxz2/2Ly 

6yc  = 

Ryc/z. 

z(m) 

k  (nr1) 

Ryc  (m) 

eyc  (rad) 

300 

420 

1.89 

0.0063 

300 

840 

3.78 

0.013 

500 

420 

5.25 

0.011 

500 

840 

10.50 

0.021 

1000 

420 

21.00 

0.021 

1000 

840 

42.00 

0.042 

We  note  that  Ryc  and  6yc  depend  on  the  parameters  {p'2},  Lx  and  Ly  through  the 
quantity  Q  =  {p'2}  Lx/Ly.  In  the  above  table,  5  =  10'7.  For  example,  if  the  anisotropy 
ratio,  Lx/Ly,  is  30  at  the  meter  scales,  rather  than  100  as  we  have  chosen  above,  all 
numbers  in  the  right-hand  two  columns  must  be  multiplied  by  .3.  We  note,  however, 
that  in  this  case  z/L,  =  2.65,  so  that  the  analytic  approximation  is  barely  adequate.  On 
the  other  hand,  if  L^Ly  =  10,  but  at  the  same  time  we  choose  {p*2}  =  10'8-  rather  than 
IQ  9,  all  numbers  in  the  table  remain  the  same  and  z/L,  =  8.82. 


F.  Point  Source  Boundary  Condition  -  Determination  of  the  Coherence  Function 

When  the  initial  plane  wave  is  replaced  by  a  point  source  (all  other  parameters 
remaining  the  same),  we  find  the  result 

1 

(r(s,z)}  =  (I jJz2)  exp(-k2z[  o(0)  -  J  o(as)  da).  (18) 

0 

Here,  lz  is  a  normalization  constant.  Apart  from  the  z2  decay,  the  above  formula  differs 

from  the  plane  wave  case  by  the  replacement  of  o(s)  by  the  integral  in  the  exponent.  It 
is  convenient  to  do  calculations  in  the  following  normalized  form: 
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(Hs.z)}  =  { r(s,z)}/(iz/z2) 


i 

=  exp[-(1  -  (1/  o(0))(  J  5(as)  da)(z/L,))],  (1 9) 

0 


where  L,  =  1/[k2  5(0)].  Here  5(0)  =  {p'2}Lx. 

Choosing  the  following  typical  values, 

k  =  420  nr1- 
Lx  =  100  m, 

{p’2}  =  10-9. 
z  =  500  m, 


we  have  z/L,  =  8.8.  This  means  we  are  in  the  strong  scattering  region.  If  Lx  and  {p'2} 
are  approximately  correct,  then  only  by  decreasing  the  signal  frequency  or  reducing 
the  range  can  we  move  out  of  this  region.  However,  for  the  previous  theory  to  be 
correct  we  require  z/Lx  to  be  at  least  greater  than  three. 

From  a  computational  point  of  view,  it  is  desirable  to  be  in  the  strong  scattering 
region  (as  we  shall  see  below).  The  disadvantage  is  that,  in  this  case,  one  can  only 
determine  the  smaller  length  scales  in  o(s).  The  larger  scales  have  little  effect  in  this 
regime. 

If  the  correlation  function  is  of  the  exponential  form 
5(as)  =  exp  (-al  sxl  /Lx  -  al  sy!  /Ly)  (20) 

we  find 

(f  (s,z)}  =  exp  {-[1  -  (1/(1  sxl  /L,  + 1  syl  /Ly  ) 

(1  -  exp-(|sxl/Lx  +  lsyl/Ly  ))](z/L,)}  (21) 


When  z/Lf  »  1 ,  only  those  values  of  s  for  which  I  sxl  /Lx  and  I  syl  /Ly  are  much  less 
than  unity  contribute  significantly  to  (f(s)}  In  this  case,  we  may  expand  the  exponential 
exp  -(I  sx!  /Lx  +  I  Syl  /Ly)  in  a  power  series  and  we  will  find,  finally 

(F(s,z)}  =  exp  (-(I  sxl  /Lx  + 1  Syl  /LyKz/L,)]. 


(22) 
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The  angular  power  spectrum  associated  with  {r(s,z)}  is  defined  as 

oo 

{ f  (k.z)}  =  //  exp  [iic  •  s]{r(s,z)}dsxdsy 

-oo 

and  we  find,  using  Eq.  (22), 

4(2Lf/z)2LxLy 

{f(K,z)}  - - .  (24) 

[(2LfKxLx/z)2  +  1]  [(2LficyLy/z)2  +  1] 

From  Eq.  (24)  we  see  that  the  characteristic  scattering  angle  in  the  vertical 
direction  is  of  the  order  of  0y/k  =  (z/LfJ^kLy).  This  result  is  the  same  as  that  obtained 
for  the  quantity  0yc  =  Ryc/z  for  the  spread  of  a  Gaussian  beam.  Thus,  we  see  that  we 
obtain  the  same  type  of  information  if  we  measure  either  the  intensity  distribution  of  a 
Gaussian  beam  or  the  coherence  function  associated  with  a  point  source.  In  view  of 
the  expected  experimental  errors,  it  would  be  most  desirable  to  measure  both 
quantities  simultaneously. 

In  the  literature,  it  is  often  postulated  that  the  wave  number  spectra  in  the 
horizontal  and  vertical  directions  obey  inverse  power  laws  of  the  form 

W(Kj)  =[r(Vj  +  V2)L/  Vic  r(Vj)] 

[1/(1  +(KjLi)2]<vi*’'2>.  (25) 

The  normalization  is  chosen  so  that  5(0)  =  1 .  Here  T(Vj)  denotes  a  Gamma  function. 
We  then  find  that  D(Sj)  =  [  5(0)  -  5(Sj)]  is  equal  to 

D(Sj)  =  [1  -  20  -V)(1/T(v,))  (s/L,)vi  Kv  (s/Lj)]  (26) 

In  the  strong  scattering  limit  we  find,  then 

(r(s,z)}  =  exp  [-(Bvx)(|  sxl  /Lx)2  +  B(vy)(|  Syl  /Ly)2  Vv)(z/Lf)]  (27) 

where 

B(Vj)  =  T(1  -  Vj)/  T(1  +  Vi)(2Vj  +  1  )2vi . 


i : 


(28) 
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For  Vj  =  V2,  we  find  Eq.  (22).  A  few  values  of  B(Vj)  are: 


Vi 

B(Vj) 

0.500 

0.500 

0.625 

0.494 

0.750 

0.558 

0.875 

0.854. 

If  Vj  =  1  then  B  diverges. 

There  is  no  simple  analytic  expression  for  the  spectrum  (r(K,z)}.  However,  it  is  a  simple 
matter  to  evaluate  the  function  numerically. 


G.  Inversion  Formulas 

The  inversion  formulas  needed  to  determine  o(s)  from  either  {H  or  {1}  are  direct 
extensions  of  the  methods  used  for  the  isotropic  case.  In  this  section,  we  shall  consider 
the  following  cases: 

1 .  Measurement  of  {T}  for  a  plane  wave  boundary  condition. 

2.  Measurement  of  {T}  for  a  point  source  boundary  condition. 

3.  Measurement  of  {1}  for  a  gaussian  beam  boundary  condition. 

4.  Measurement  of  {R|}  for  a  plane  wave  boundary  condition. 

1.  Plane  wave  boundary  condition  for  (H: 

Direct  solution: 

{T(s,z)}  =  I  exp(-k2[o(0)  -  c(s)]  z)  (29) 

Inverse  solution: 

o(s)  =  5(0)  +  (1/k2z)  log[{T(s,z)}/I]  (30) 

or 


1  3 


o(s)  =  (1/k2z)  log({r(s,z)}/{r(oo,z)}]. 


(31) 
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2.  Point  source  boundary  condition  for  (D: 

Direct  solution: 

1 

{r{s,z)}  =  (l2/z2)  exp(-k2z[5(0)  -  Ja(as)da])  (32) 

0 

Inverse  solution: 

1 

J  5(as)da  =  5(0)  +  (l/k2z)  logKHs.zM^z2)].  (33) 

0 

The  coordinates  sx  and  Sy  may  be  any  orthogonal  coordinates  in  the  transverse 

plane.  To  reduce  Eq.  (33)  to  an  expression  for  5(s),  we  choose  one  of  the  s 
coordinates  equal  to  zero.  For  example,  we  take  here  sy  =  0.  Next,  noting  that 

1 

d/dsx[  sx  j  5(asx)da]  =  o(sx) 

0 

we  find  the  following  expression  for  5(sx) 

5(sx)  =  5(0)  + 

d/dsx[(sx/k2z)  log({r(sx,z)}/(lz/z2))]  (34) 

or 

5(sx)  =  d/dsx[(sx/k2z)  log({r(sx,z)}/{r(co,z)})].  (35) 

There  is  a  similar  expression  for  sy  when  sx  =  0. 

3.  Gaussian  source  boundary  condition  for  (11: 

Direct  solution: 

oo 

(l(p,z)}  =  (I  b2/4ji)  JJ  dw  exp(-iw  •  p)  exp[-w2b2(1  +  z2/k2b4)/4  +  Q(w,z)]  (36) 


T 
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I 


where 


wz/k 


Q(w,z)  =  (k3/w)  J  [  5(wxt/w,wyt/w)  -  5(0,0)]dt. 


(37) 


Here,  {l(p,0)}  =  I  exp[-p2/b2]. 

Inverse  solution: 


Q(w,z)  =  w2b2{1  +  z2/k2b4)/4  +  log[1/(l  b2x)  JJ  {l(p,z)}  exp(iw  •  p)dp],  (38) 


To  find  5  from  Q  along  one  of  the  transverse  axes,  we  set  one  of  the  orthogonal 
coordinates  of  w  equal  to  zero,  say  wy  =  0.  We  find,  then,  from  Eq.  (37) 

5(wxz/k,0)  =  5(0,0)  +  (1/k2z)[5(wxQ(wx,0,z))/3wx].  (39) 

There  is  a  similar  expression  for  wy  when  wx  =  0.  Moreover,  by  suitable  rotation 
of  the  transverse  coordinate  axes,  we  may  find  5  along  any  axis.  The  coupling 

between  coordinate  axes  in  5  may  be  obtained  by  keeping  w  arbitrary  and 
differentiating  Eq.  (37)  with  respect  to  the  vector  magnitude,  w.  The  expression  is  then 


I 

I 

I 

I 

I 


5(  wz/k)  =  5(0,0)  +  (1/k2z)[9(wQ(w,z))/dw]. 

4.  Plane  wave  boundary  condition  for  (FM: 

Direct  Solution: 

The  perturbation  solution  for  {R|(xT1,xT2,z)}  is 

oo 

(R|(XT1,XT2,Z)  =  l2[1  +  (k2z/2)  Jj  5(k12)  exp(ik12  •  r12) 

-oo 

(1  -  (k/k122z)  sin[k122z/kj)dk12], 


I  s 


(40) 


(41) 
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wh  3re 


oo 

o(k12)  =  (1/71)2  JJ  a(r12)  exp(ik12  •  r12)dr12 

-oo 


From  Eq.  (41).  we  see  that  (R|)  is  determined  by  a  rather  complicated  integral  of 
the  function  S(k12).  However,  this  integral  reduces  to  a  simplified  relation  in  two  limits 
(1)  in  the  geometric  optics  region  where  z  «  ksc2  (sc  is  the  smallest  characteristic 
length  associated  with  the  coherence  function)  and,  (2)  when  z  »  ksc2.  In  the 
geometric  scattering  region  we  find 


oo 

{R|(Xti,xT2.z)  =  I2  [1  +  (V12)(z3)  JJ  o(k12)(k124)  exp(ik12  •  r12)dk12], 

-oo 

=  |2[i  +(V3)z3y45(r12)]. 

When  z  »  ksc2, 

oo 

(R,(s,z)}  =  I2  [1  +  (k2z/2)  JJ  5(k12)  exp(ik12  •  r12)dk12] 

-oo 

=  l2[1  +  2k2z  o(r12)]. 

Inverse  solution: 

For  z  «  ksc2  we  find 

V4  5(r12)  =  (3/z3)[Rl(r12,z)  -  l2]/l2. 

For  z  »  ksc2  ,we  have 
5(r12)  =  (1/2k2z)(R,(r12,z)  - 12]/!2. 


(42) 


(43) 


(44) 


(45) 
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III.  Application  of  the  Theory  to  the  Arctic  Environment 


To  properly  design  and  evaluate  an  experiment  to  be  conducted  under  Arctic  ice  , 
it  is  desirable  to  have  some  idea  of  the  scattering  properties  of  the  Arctic  Ocean 
environment.  Despite  an  extensive  search  of  the  literature,  very  little  data  was  found 
for  the  meter-to-centimeter  scales  in  which  we  are  interested.  Not  only  is  the  shape  of 
the  correlation  function  (or  associated  wave  number  spectrum)  largely  unknown,  but 
the  relevant  length  and  intensity  scales  can  be  estimated  at  best  only  to  within  an 
order  of  magnitude.  The  best  information  that  we  have  been  able  to  obtain  is  given  in 
the  references  at  the  conclusion  of  this  section. 

For  purposes  of  estimating  effects  of  the  fine  structure  in  the  ocean  we  have 
assumed  the  values  of  the  ratio  {c'2}/{c}2  (the  ratio  of  mean  square  sound-speed 
flucuations  to  the  mean  sound-speed  squared)  to  be  of  the  order  of  10-10  to  10*9.  The 
vertical  correlation  length  is  taken  to  be  1-3  meters,  with  the  horizontal  correlation 
length  about  10  to  100  times  the  vertical  correlation  length.  Typical  results  using  these 
numbers  are  given  in  Section  II. 
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IV.  Effect  of  Random  Velocity  Fluctuations  on  Underwater  Scattering 


There  is  a  well-developed  theory  to  account  for  the  effect  of  random  fluctuations 
in  the  sound-speed  on  the  underwater  scattering  of  acoustic  radiation.  In  Appendix  A, 
we  show  how  the  theory  may  be  extended  to  include  the  effect  of  random  fluctuations 
in  fluid  velocity.  We  assume  that  the  characteristic  time  and  length  scales  associated 
with  the  velocity  correlation  functions  are  long  when  compared  to  the  corresponding 
scales  (period  and  wave  length)  of  the  acoustic  radiation. 

The  importance  of  fluid  velocity  fluctuations  depends  upon  the  scale  sizes  of 
importance  in  any  acoustic  experiment.  For  experiments  in  which  internal  waves  have 
the  dominant  effect  on  acoustic  fluctuations,  the  fluid  velocity  fluctuations  may  be 
neglected.  However,  for  experiments  in  which  the  fine  structure  plays  a  dominant  role, 
fluid  velocity  fluctuations  begin  to  have  an  effect.  In  the  paper  included  as  Appendix  A 
we  show  how  the  fluid  velocity  fluctuations  may  be  taken  into  account  using  the 
inverse  method  proposed  for  the  arctic  experiment.  More  standard  methods  require 
direct  measurement  of  the  fluid  velocity  fluctuations;  a  procedure  that  is  exceptionally 
difficult  and  rarely  reported  in  the  literature. 


V.  Determination  of  the  Intensity  Distribution  from  a  Circular  Source  in  an  Isotropic 
Random  Medium 

In  a  paper  by  Whitman  and  Beran  (1970),  the  intensity  distribution  was 
determined  for  a  Gaussian  beam  propagating  in  an  isotropic  random  medium.  Here, 
we  present  a  similar  result  for  a  circular  source  (i.e.  constant  intensity  for  a  radial 
distance  less  than  the  radius  a  and  zero  for  greater  radial  distances).  The  difficulty  in 
the  calculation  is  determining  the  initial  coherence  function,  since  the  initial  field  is 
discontinuous.  We  shall  first  outline  the  method  of  solution  and  then  determine  the 
initial  coherence  function.  Finally,  we  present  the  integral  that  must  be  evaluated  to 
determine  the  intensity. 


A.  Method  of  Solution 

The  basic  equation  for  the  coherence  function,  (T),  is 

a{r}/az  =  (i/k)  a/as  •  ap(n  +  k2F<s){n,  (46) 

where  z  is  the  range,  p  is  the  mean  transverse  coordina'e,  s  is  the  transverse 
difference  coordinate  and  k  is  the  wave  number.  The  function  F(s)  is  given  by  the 
expression 

F(s)  =  [o(s)  -  5(0)] 

where  o(s)  is  the  transverse  correlation  function  associated  with  the  index-of-refraction 
fluctuations. 
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To  solve  Eq.  (46),  we  take  the  double  Fourier  transform  with  respect  to  p.  The  result  is 

3{r}/az  -  (i/k)u  •  a/as(n  ■  k2F(s)(n  =  o.  (47) 

We  solve  Eq.  (47)  by  the  method  of  characteristics.  The  result  is,  (after  some 
manipulation), 

z 

(l(u,z)}  =  (r(uo,so,0)}  exp[k2  J  F(uq/k)dq]  (48) 

0 

where  s0  =  -uGz/k  and  u  =  uQ.  The  intensity  is  found  from  (l(u,z)}  from  the  Bessel 
transform  relation 

OQ 

{l(R.z)}  =  (1/271)  J  (I(U.Z)}  J0(uR)udu.  (49) 

0 

B.  Determination  of  the  Initial  Coherence  Function 
In  terms  of  the  initial  field,  V(x),  we  have 

oo 

(r(uo,so,0)}  =  jj  exp[iu„  •  p0]  V(pQ  +  s0)  V'(p0)dp0  (50) 

-oo 

Using  the  convolution  theorem,  this  yields 

oo 

{T(uo,so,0)}  =  (V2n  )2  jj  exp[-s0  •  t]W(t)  W*(t  +  u0)dt,  (51 ) 

-oo 

where 

oo 

W(t)  =  jj  exp[it  •  p0]  V(p0)dp0  (52) 

-oo 

For  a  circular  source,  we  find 

W(t)  =  lQ  2rea  J, (at )/t.  (53) 
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Since  the  function  {r(uo,so,0)}  depends  only  on  u0,  we  use  the  simplified  notation 
ls(u0,o)  =  {r(uo,so,0)} 
and  find,  finally, 

°°  2tc 

ls(u0,o)  =  a2  JJ  exp[iu0t  cos(v)  z/k]  J^atyt 
00 

J^a^2  +  u02  +  2t  u0cos(v))1/2] 

1/(t2  +  u02  +  2t  u0  cos(v))1/2dt  dv.  (54) 

The  final  solution  is  then 

z 

(l(u,z)}  =  ls(u,0)  exp[k2  |  F(uq/k)dq],  (55) 

0 

where  ls(uo,0)  is  given  by  Eq.  (54)  (here  u  =  u0)  and  (l(R,z)}  is  then  found  using 
Eq.  (49). 


We  note  that,  in  the  far  field  z  -» °°,  only  values  of  Uo/k«1  are  important,  and  we 
determine  that 


ls(u0,o)  =  l0  27ta2  J  [J-|(at)/t]2  J0(u0tz/ak)tdt. 
0 


(56) 
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VI.  Studies  of  the  Optimum  Parameters  for  a  Tank  Experiment 


As  a  part  of  our  work,  we  attemped  to  conduct  a  limited  validation  of  the  theory 
presented.  In  order  to  acomplish  this  within  the  scope  of  this  current  contract,  we 
performed  an  experiment  in  which  we  scattered  sound  from  artificially  generated 
sound  speed  inhomogenieties  in  a  large  laboratory  tank.  In  designing  this  experiment, 
we  needed  to  determine  the  correlation  function  associated  with  the  sound-speed 
fluctuations,  using  the  inverse  method.  In  order  to  do  this,  the  effect  of  scattering  must 
be  stronger  than  source  diffraction  effects,  but  not  so  strong  that  only  the  very  smallest 
scales  of  the  fluctuations  are  important.  The  optimum  conditions  occur  when  the  the 
strength  parameter  B  defined  below,  is  of  order  unity  and  the  ratio  s/a  (scattering  to 
diffraction  effect)  is  as  large  as  possible. 


A.  Temperature  Fluctuations 

While  the  ultimate  parameter  of  interest  was  the  sound  speed,  cost 
considerations  led  us  to  measure  the  temperature  of  the  water  and  then  to  compute 
sound  speed  for  zero  salinity  (tap)  water.  The  relationship  between  the  temperature 
and  sound  speed  fluctuations  may  be  determined  from  the  following  well  known 
relation  between  sound  speed  c  and  temperature  T. 

c  =  1449.0+  4.6  T  -  0.055  T2  +  0.0003  T3, 

where  c  is  in  m/sec  and  T  is  in  degrees  Centigrade.  Taking  the  differential  of  the 
above  equation  yields 

dc  =  4.6  dT  -  0.1 1  dT  +  0.0009  T2  dT. 

For  T  =  25°  C,  a  typical  tank  temperature, 


dc  =  2.25  dT. 

For  dT  =  1°  C,  dc/c  =  0.0015  and  p2  =  4(dc/c)2  =  9x10-6=  10-5;  p2  =  {^2}. 

B.  Strength  of  Scattering 

The  strength  of  scattering  B  may  be  defined  as  follows 
B  =  p2k2/z 

where  k  is  the  radiation  wavenumber  (m-1>and  for  a  frequency  of  f=1  MHz, 

(k  =  4200  nr1) ,  /  is  the  characteristic  turbulent  correlation  length  and  z  is  the  path 
length  (m). 
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Typical  values  are  as  follow. 


k(mA) 

F2 

/(m) 

z(m)  B 

B1/2 

4200 

10-5 

.1 

3 

52.90 

7.27 

4200 

10-5 

.01 

3 

5.29 

2.30 

2100 

10-5 

.01 

3 

1.32 

1.15 

2100 

10-5 

.03 

3 

3.97 

1.99 

4200 

10-6 

.03 

3 

1.59 

1.26 

2100 

IQ'6 

.03 

3 

0.40 

0.63 

C.  Angle  of  Scattering  and  Spread 

Single  scattering  angle,  0TS  =  (1/k/)- 


Multiple  scattering  angle,  0T  =  B1/2(1/k /). 


Spread,  s  =  (0T)z. 


D.  Fresnel  Region 

For  a  Gaussian  source,  the  diffraction  effect  is  a  minimum  at  a  distance  z  if  the 
source  size  a  is  chosen  so  that  a  =  (z/k)1/2  . 

Typical  values  are  shown  below. 

z(m)  kfrrr1)  a(m)  ka 


3 

3 


4200 

2100 


0.0267 

0.0378 


112.2 

79.4 
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E.  Ratio  s/a 


~r 


The  ratio  s/a  =  [Bz/k/2]1'2  gives  us  the  relative  importance  of  turbulent  scattering  to 
source  diffraction.  Some  typical  values  are  given  below. 


Table  I 


kfnr1) 

H2 

l(  m) 

z(m) 

B 

Q^/2 

s/a 

8400 

10-5 

0.1 

3 

211.600 

14.540 

2.740 

8400 

10'6 

0.1 

3 

21.600 

4.600 

0.866 

8400 

10-7 

0.1 

3 

2.120 

1.450 

0.274 

8400 

10-5 

0.03 

3 

63.600 

8.000 

7.140 

8400 

10-6 

0.03 

3 

6.360 

2.520 

2.250 

8400 

10-7 

0.03 

3 

0.636 

0.800 

0.710 

8400 

10-5 

0.01 

3 

21.160 

4.600 

8.680 

8400 

10-6 

0.01 

3 

2.120 

1.450 

2.740 

8400 

10-7 

0.01 

3 

0.212 

0.460 

0.866 

4200 

10-5 

0.1 

3 

52.900 

7.270 

1.940 

4200 

10-6 

0.1 

3 

5.290 

2.300 

0.613 

4200 

10-7 

0.1 

3 

0.529 

0.727 

0.194 

4200 

10  s 

0.03 

3 

15.900 

3.990 

5.020 

4200 

10-6 

0.03 

3 

1.590 

1.260 

1.590 

4200 

10-7 

0.03 

3 

0.159 

0.399 

0.502 

4200 

10-5 

0.01 

3 

5.290 

2.300 

6.140 

4200 

10-6 

0.01 

3 

0.529 

0.727 

1.940 

4200 

10-7 

0.01 

3 

0.053 

0.230 

0.614 

2100 

10-5 

0.1 

3 

13.200 

3.630 

1.370 

2100 

10-6 

0.1 

3 

1.320 

0.363 

0.137 

2100 

10-5 

0.03 

3 

3.970 

1.990 

2.510 

2100 

10-6 

0.03 

3 

0.400 

0.630 

0.800 

2100 

10-7 

0.03 

3 

0.040 

0.200 

0.253 

2100 

10-5 

0.01 

3 

1.320 

1.150 

4.340 

2100 

10-6 

0.01 

3 

0.132 

0.363 

1.370 

2100 

10-7 

0.01 

3 

0.013 

0.115 

0.434 

The  best  situation,  from  an  inversion  point  of  view,  is  B  near  unity  and  s/a  as  high 
as  possible.  As  we  stated  above,  the  reason  to  require  that  the  value  of  B  is  close  to 
one  is  to  insure  that  all  scales  of  the  turbulence  contribute  to  the  intensity  pattern.  If  B 
is  very  large,  only  the  very  small  scales  are  important.  On  the  other  hand,  if  B«1 ,  then 
the  turbulence  does  not  have  much  effect. 
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VII.  Laboratory  Tank  Experiment 

A  prototype  experiment  was  conceived  and  executed  with  the  objective  of  at  least 
partially  validating  the  theory  presented  in  the  preceding  part  of  this  report.  The 
experimental  protocol  was  heavily  influenced  by  the  availability  of  certain 
instrumentation  and  computer  data  acquisition  interfaces.  Only  with  the  use  of  those 
resources  could  we  hope  to  accomplish  the  objective  of  comparing  theory  and 
experiment  within  our  existing  funding  resources. 

Basically,  the  experiment  was  designed  in  two  parts.  The  first  involved 
measurement  of  the  correlation  function  for  an  artificially  induced  heterogeneous 
spatial  structure  in  the  speed  of  sound.  The  structure  was  created  in  two  ways.  One 
involved  the  placement  of  a  heater  grid  on  the  tank  bottom.  A  thin  plastic  sheet  was 
placed  over  this  heater  grid.  Holes  had  been  punched  (ca  1  cm  diameter)  randomly  in 
the  sheet  in  a  position  that  was  directly  under  the  projected  sound  propagation  path. 
The  density  of  these  holes,  through  which  the  heated  water  rose  (due  to  its  own 
density  contrast  with  the  surrounding  water)  was  2000  /  m2.  The  second  technique 
involved  a  reversal  of  the  process  described  above.  Ice  was  floated  in  the  water 
above  a  second  sheet  that  was  similiar  to  that  described  for  the  heater.  As  the  cold 
water  sank,  it  was  introduced  into  the  propagation  path.  Not  suprisingly,  as  the 
experimental  tank  is  outdoors,  the  heater  resulted  in  the  strongest  sound  speed 
heterogeneity  in  the  experiments  conducted  in  the  winter  when  the  ambient  tank 
temperature  was  about  8°C,  and  the  melting  ice  was  most  effective  in  the  summer  and 
fall  when  ambient  air  temperatures  were  higher  (ca  20-23°C).  On  occasion,  when  the 
ambient  temperature  was  intermediate  (ca  14°C),  and  a  strong  scattering  field  was 
desired,  both  techniques  were  used  simultaneously. 

A.  Sound  Speed 

The  sound  speed  was  computed  from  measurements  of  the  water  temperature  at 
six  points  separated  by  1  cm  along  a  horizontal  line.  Since  tap  water  was  used  in  the 
tank,  the  salinity  was  assumed  to  be  sufficiently  close  to  zero  to  justify  an  assumption 
that  it  could  be  neglected  both  on  the  average  and  in  terms  of  the  spatial  heterogeneity 
that  it  might  introduce.  The  thermistors  used  were  calibrated  to  an  accuracy  of  0.01  °C 
in  an  oil  bath.  The  time  constant  for  these  thermistors  had  been  previously  measured 
and  were  found  to  be  ca  8  ms  in  water  (1/e).  A  regulated  voltage  source  was  used  to 
excite  each  thermistor  and  the  resulting  voltages  were  digitized  by  independent  auto¬ 
zeroing  A/D  converters.  The  circuit  was  designed  to  automatically  compensate  for  any 
fluctuations  in  the  source  voltage  by  use  of  a  ratiometric  technique. 

The  A/D  converters  (12-bit)  were  strobed  to  begin  their  conversions 
simultaneously  and  the  digitized  data  were  latched  for  transfer  to  the  data  acquisition 
computer  as  soon  as  each  set  of  conversions  was  completed.  Thus,  a  complete  set  of 
independent  temperatures  was  available  at  a  30  Hz  rate.  These  binary 
representations  of  the  temperature  at  each  thermistor  were  transfered  to  the  data 
acquisition  computer  as  a  serial  sequence  of  bytes.  There  they  were  stored  in  a  RAM 
buffer  until  the  required  data  sequence  was  completed.  At  that  time  the  numbers  were 
converted  to  engineering  units  (°C)  and  stored  on  a  floppy  disk.  An  option  was 
provided  for  deleting  unwanted  data  from  the  data  stream  before  storing  in  the  buffer  if 
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the  full  sampling  rate  of  30  Hz  was  not  desired.  This  procedure  was  implemented  by 
the  simple  process  of  skipping  the  storage  of  data  sets  between  desired  sample  times. 

The  data  were  displayed  in  both  a  columnar  format  on  a  CRT  and  optionally  on  a 
printer  and  as  a  plot  on  a  graphic  display.  The  temperature  data  were  then  converted 
to  sound  speed  and  the  correlation  function  for  that  parameter  was  also  made 
available  in  both  numerical  and  graphical  display  formats. 

B.  Fluctuations  in  the  Received  Acoustic  Pulses 

Acoustic  pulses  were  projected  over  a  path  of  2.33  m  at  a  depth  of  1  m  in  the  tank. 
The  water  depth  was  maintained  at  just  under  2  m.  The  tank  is  2  m  (deep)  by  2  m 
(width)  by  4  m  (length).  A  pulse  length  of  1  ms  was  used  and  the  acoustic  center 
frequency  was  0.925  MHz.  The  projector,  a  transducer  constructed  from  a  PVDF 
(polyvinylidene  fluoride)  film  of  1 10  4m  thickness  was  fabricated  as  five  concentric 
rings.  The  outer  diameter  of  the  largest  ring  was  10  cm.  The  drive  voltage  to  these 
rings  was  shaded  in  an  approximation  to  a  Gaussian  curve  with  a  resistive  network. 
The  resultant  beam  pattern  was  a  good  approximation  to  a  Gaussian  shape  inside  the 
Ve  points  on  the  intensity  pattern.  Beyond  that  point,  however,  the  approximation 
deviated  from  a  Gaussian  shape,  falling  off  too  rapidly.  It  is  possible  that  this  deviation 
was  due  to  an  uncompensated  reactive  component  of  the  cable  and  shading  network 
at  this  (relatively  high)  frequency.  This  kind  of  error  is  layout  sensitive,  and  can  not  be 
modeled  to  allow  ready  correction.  It  can  be  adjusted  by  an  iterative  approach, 
however,  our  time  and  resources  did  not  permit  this  to  be  accomplished  for  this  task. 

Two  small  probe  hydrophones,  one  a  disk  (ca  2.5  cm  diameter  and  the  other  a 
cylinder  (1.1  cm  height  and  0.7  cm  outer  diameter)  were  located  at  the  same  depth  as 
the  source  at  a  distance  of  2.33  m.  The  disk  was  fixed  relative  to  the  transmitted  sound 
beam  and  the  cylinder  could  be  adjusted  in  the  horizontal  in  a  direction  that  was 
orthogonal  to  the  undistorted  direction  of  travel  of  the  source’s  main  lobe. 

In  any  acoustics  experiment  to  be  conducted  in  a  finite  sized  tank,  one  must  be 
concerned  with  the  potential  impact  of  reflections  on  the  received  signals.  In 
particular,  at  the  frequencies  with  which  we  were  concerned  in  this  experiment, 
modulation  of  the  surface  reflected  path  by  small  ripples  would  have  interfered  with 
our  ability  to  discriminate  forward  scattering  induced  reflections.  In  order  to  assure  that 
surface  reflections  would  not  be  present  in  our  measurements,  we:  1)  chose  to  work 
with  as  short  a  pulse  as  was  practical  (<  1  ms),  2)  used  a  range  gate  to  assure  that  the 
measurement  included  only  the  direct  arrival  and  those  (relatively)  small  variations 
due  to  forward  volume  scattering  rather  than  signals  reflected  from  the  water  surface, 
the  tank  walls  or  bottom,  3)  used  the  directivity  of  the  receiving  hydrophones  to  reject 
energy  arriving  at  the  range  gate  during  the  sampling  period  at  large  vertical  angles  off 
their  pointing  direction  and,  4)  transmitted  at  a  rate  (ca  1/sec)  that  assured  that  the 
reverberation  in  the  tank  from  previous  transmissions  had  been  sufficiently  attenuated 
before  the  next  sequence  began. 

A  bandpass  filter  centered  at  0.925  MHz  was  used  to  reject  electronic  noise 
(principally  from  commercial  radio  stations  in  this  band)  and  thermal  noise  that  fell 
outside  the  desired  band  around  the  acoustic  transmission  frequency.  The  energy  in 
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VIII.  Analysis  of  the  Tank  Experiment  Results 


The  experimental  results  obtained  demonstrate  that  under  laboratory  conditions, 
it  is  feasible  to  measure,  both  the  correlation  function  of  the  sound-speed  fluctuations 
and  the  intensity  distribution  of  a  beam  that  has  traversed  a  heterogeneous  medium. 
This  means  that  it  is  possible  both  to  study  the  validity  of  the  basic  theory  that  has  been 
developed  to  predict  the  intensity  distribution  and  to  determine  the  utility  of  the  inverse 
procedure  discussed  in  the  previous  sections.  We  have  shown  that  it  is  possible  to 
obtain  values  of  p2  of  the  order  of  10'7  in  the  laboratory.  For  this  value  of  p2,  significant 
scattering  of  a  beam  occurs  and  the  scattering  effect  may  be  distinguished  from  the 
diffraction  effect.  In  addition,  it  has  been  possible  to  measure  the  correlation  of 
intensities  so  that  measured  values  of  the  scintillation  index  may  also  be  compared  to 
those  predicted  by  theory. 

An  example  of  a  measured  sound-speed  correlation  function  is  given  in  Fig.  3. 
The  characteristic  correlation  length  is  about  1  cm  and  p2  =  0.806x1 O'6.  By  changing 
the  heating  or  cooling  in  the  tank  p2  may  vary  between  10  ®  and  106.  Thus,  we  can 
meet  the  conditions  that  B  is  approximately  unity  and  s/a  >  1 .  See,  for  example, 
k  =  4200  nr1 ,  /  =  0.01  and  0.03,  p2  =  1 0'6  in  Table  I. 


d(cm) 

Fig.  3:  Sound  Speed  Correlation  Function  vs.  Separation  Distance 


In  Table  II,  we  show  normalized  intensity  distributions  that  were  measured  under 
weak  turbulence  conditions  (p2  =  0.5x1  O'7)  and  strong  turbulence  conditions 
(  p2  =  0.806x1  O'6).  We  see,  from  comparison  of  the  intensity  distributions,  that  the 
beam  in  the  strong  turbulence  case  was  wider.  In  this  feasibility  experiment,  the 
absolute  intensities  were  not  measured  (though  this  could  be  done  by  calibrating  the 
probe  hydrophones  by  a  reciprocity  or  a  comparison  technique).  We  did  not  measure 
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the  intensity  reduction  on  the  beam  axis  in  the  current  effort.  The  beam  shape 
departed  too  far  from  Gaussian  beyond  the  1/e  -  point  to  make  it  useful  to  carry  out  an 

inverse  calculation  on  these  data. 

In  future  work,  we  have  two  choices.  We  can  either  devote  considerable  more 
effort  to  making  the  skirts  of  the  beam  distribution  more  nearly  Gaussian  (a  tedious,  but 
achievable  goal)  or  we  can  fully  develop  a  theory  for  a  beam  with  a  more  general 
beam  shape.  Since  the  natural  shape  of  lines  and  disks  involve  the  appropriate 
Bessel  functions,  a  model  for  such  a  beam  was  considered  on  a  theoretical  basis  and 
seems  practical.  However,  the  numerical  calculations  tor  a  non-uaussian  beamwas 
not  within  the  scope  of  the  current  work.  In  either  case,  the  theory  may  then  be  directly 
compared  to  experiment  and  we  may  further  test  the  inverse  procedure.  Here^  we  note 
only  that  the  measured  intensity  distribution  gives  results  that  are  of  the  correct  order  o 
magnitude.  For  example,  the  mean  free  path  of  scattering,  lF,  may  be  defined  for  this 
experiment  as  lF  =  (^k2/)1.  For  strong  scattering,  we  find  lF  is  approximately  6  m.  (In 
the  experiment  k  =  4000  nr’).  The  propagation  path  was  2.33  m  so  that  we  expect  that 
roughly  30%  of  the  beam  was  scattered.  In  addition,  the  transverse  scattered  distance 
at  2.33  m  was  (2.33)/(k/>  ,  which  is  approximately  6  cm.  These  numbers,  which  are 
order  of  magnitude  estimates,  are  consistent  with  the  results  of  Table  II. 


Table  II 


Distance  Intensity  Intensity 

from  beam  axis  Weak  Scattering  Strong  Scattering 
(cm)  (Normalized)  (Normalized) 

0  1.000  1  00 

1  0.910  0.96 

2  0.700  0.87 

3  0.360  0.63 

4  0.100  0.28 

6  0.003+  0.034 

8  0.011  0.032 


+  near  measurement  system  noise  limit. 
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IX.  Conclusions  and  Recommendations 

The  basic  theory  required  for  the  development  of  a  new  acoustic  technique  for 
studying  fine-to-intermediate  scale  processes  and  distributions  of  properties  in  the 
ocean  have  been  successfully  addressed.  Insofar  as  we  were  able  to  test  that  theory, 
experiment  and  theory  are  in  agreement.  A  critical  need  exists  to  understand  ocean 
processes  on  scales  of  100's  of  meters  to  centimeters.  The  tools  and  techniques  that 
are  currently  available  are  inadequate.  There  are  potential  applications  of  such  a  tool 
in  addressing  key  problems  in  underwater  acoustics,  and  in  physical,  chemical  and 
biological  oceanography.  We  conclude,  based  on  our  results  to  date,  that  additional 
work  is  justified  in  pursuing  this  approach  to  a  better  description  and  understanding  of 
intermediate-to-fine  scale  oceanography. 
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APPENDIX  A 


Effect  of  Random  Velocity  Fluctuations 
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Abstract 

There  is  a  well-developed  theory  to  account  for  the  effect  of  random  fluctuations  in  the 
sound-speed  on  the  underwater  scattering  of  acoustic  radiation.  Here  we  show  how 
the  theory  may  be  extended  to  include  the  effect  of  random  fluctuations  in  fluid  velocity. 
We  assume  that  the  characteristic  time  and  length  scales  associated  with  the  velocity 
correlation  functions  are  long  when  compared  to  the  corresponding  scales  (period  and 
wavelength)  of  the  acoustic  radiation.  We  conclude  with  a  discussion  of  the  application 
of  the  theory  to  ocean  studies  and  the  utility  of  using  an  inverse  method  to  find  the 
lateral  correlation  function. 


Introduction 

In  two  recent  articles1-2,  parabolic  equations  were  derived  to  take  into  account  the 
effect  of  ocean  currents  on  sound  transmission.  Attention,  in  these  papers,  was  limited 
to  the  effect  of  horizontal  ocean  currents.  However,  their  formulation  is  easily  extended 
to  include  a  full  three-dimensional  velocity  vector  field,  v(x).  Here  we  are  interested  in 
the  equations  governing  the  coherence  functions  of  the  acoustic  radiation  and  shall 
restrict  our  attention  to  random  velocity  fields  that  satisfy  the  following  conditions: 

a)  |v|/c0  «  1 

b) coTc  »  1 

c)  kL,. »  1 

where  c0  is  a  characteristic  speed  of  sound,  o>  is  the  radian  frequency  of  the  acoustic 
wave,  k  is  the  radiation  wavenumber,  Tc  is  the  shortest  characteristic  time  associated 
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with  the  velocity  correlation  functions  and  !_<.  is  the  smallest  characteristic  length 
associated  with  the  velocity  correlation  functions. 

With  the  above  restrictions,  it  is  possible  to  show  that  for  a  cw  signal  the  spatial 
part  of  the  pressure,  p(x),  satisfies  the  equation 

V2p  +k2p(x)  +  2ikoV(x)  *  Vp(x)/c  =  0.  (1 ) 

In  the  derivation  of  Eq.  (1)  we  used  the  conditions  (a)-(c).  It  was  not  necessary  to 
assume  that  the  magnitude  of  the  vertical  component  of  the  velocity  was  much  smaller 
than  that  of  the  horizontal  components. 

To  derive  a  parabolic  equation  we  write 

P(x,y,z)  =  u(x,y,z)  expflkoZ), 

where  z  is  the  mean  propagation  direction,  y  the  depth  direction,  x  the  cross-range 
direction  and  ko  =  oo/c0.  The  first  two  terms  of  Eq.  (1)  yield  the  usual  parabolic 
equation.  The  last  term  gives 


2ikQ[vT  *VTp  +  vzdu/3z  (exp(ikz))  +  ikoVzp]/c0 


In  the  appendix,  we  show  that  only  the  last  term  contributes  to  dominant  order  in 
the  parabolic  equation.  We  thus  find 

2iko3u/5z  +  V^u  +  [k2  -  k*,2  -  2k02vz/c0]u  =  0.  (2) 

When  the  sound-speed  fluctuations  are  small  we  may  write 
k2- 1^2  =  ^[2(0  -c0)/co  +  2c7co], 
where  c  is  the  mean  speed  of  sound.  This  yields 

2ikQdu/dz  +  V-^u  +  ko2[  p(x)  +  p'(x)]u  =  0  (3) 

where  p(x)  =  2(c  -  c0)/c0  and  p’(x)  =  2(c'  -  vz)/c0. 

Eq.  (3)  is  exactly  the  same  equation  as  that  used  in  studies  of  scattering  of  acoustic 
energy  by  sound-speed  fluctations  alone.  The  only  difference  is  that  here  the  velocity 
component  in  the  range  direction,  vz,  has  been  added  to  c\  Therefore  the  entire 
statistical  theory  used  to  study  scattering  from  sound-speed  fluctuations  alone  may 
also  be  used  when  we  wish  to  consider  the  effect  of  a  random  velocity  field.  In  the  next 
section  we  shall  briefly  review  the  relevant  statistical  equations  and  in  the  final  section 
we  shall  discuss  the  application  of  the  theory  to  real  ocean  studies. 
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Statistical  Theory 

To  treat  the  problem  of  scattering  from  a  stochastic  point  of  view,  we  define  the 
following  coherence  functions 


{r(xT1  ,xT2,z)}  =  {u{xT1  ,z)u‘(xT2,z)}  (4) 

{L1(xT1,xT2>xT3)xT4,z)}  =  {u(xT1,z)u*(xT2lz)u*(xT3,z)u(xT4lz)}.  (5) 

Here  u*  denotes  the  complex  conjugate  of  u  and  the  vector  xT  has  the  transverse 
coordinates  x,y.  The  brackets,  {  },  denote  an  ensemble  average. 

The  intensity,  {l(xT,z)},  is  found  from  {r(xT1,xT2,z)}  by  setting  xTl  =  xT2.  The  intensity 
fluctuations,  (l2(xT,z)},  and  the  coherence  of  intensities,  {R|(xT1,xT2,z)},  are  found  from 
{L1}  from  the  following  relations 


{I2(xt,z)}  =  {L1  (xT,xT,xT,xT,z)  (6) 

{R|(XT1.X72.2)}  -  {L1(Xti.Xti.XT2,X72.Z)}.  (7) 


We  write  p(x)  in  the  form 

M-(x)  =  p(x)  +  p'(x),  (8) 

where  jl(x)  represents  the  mean  background  profile  and  p'(x)  represents  the  random 
fluctuation  field.  We  shall  assume  here  that  the  statistics  of  the  field  are 
homogeneous.  The  two-point  correlation  function,  o(x1,x2),  is  defined  as 

a(x1,x2)  =  {p’(x1)p,(x2)}  (9) 

and  depends  upon  r  =  -  x2.  In  addition  we  define  the  function 


oo 

5(xTi  -  xT2)  =  (V4)  j  a(xTi  -  xT2,rz)drz 

-  oo 


(10) 


The  equations  governing  {r}  and  {L1}  (See  Refs.  3  -  7)  are 
3{r}/9z  =  (i/2ko)[V2Ti  -  - 


ko2(o(0)  -  5(xT1  -  xT2)]{r>  +ik0[p(xT1)  -  p(xT2)](r} 


(11) 
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a{U}/9z  =  (i^koKyZr,  -V2T2  -V2T3  +72t4]{L1}  + 


ko2[o(xTi  -  xT2)  +  o(xT1  -  xT3)  + 


o(Xf2  -  XT4)  +  5(x-p3  -  XT4)  - 


5{Xji  -  xT4)  -  a(XT2  -  x^)  - 


2a(0)]{L1}  +  ikoP(xTi)  - 

P(XT2)  ‘  P(XT3)  +  P^T^K^-1}  (1  2) 


In  the  literature,  there  are  many  solutions  of  Eqs.  (11)  and  (12)  under  a  variety  of 
conditions  (see,  for  example,  Refs.  3-18).  It  is  not  our  purpose  to  review  these 
solutions  here,  but  rather  to  emphasize  that  any  solution  now  available  can  be  used  to 
include  the  effects  of  random  velocity  fluctuations  provided  that  the  correlation  function 
a(x1tx2)  is  defined  as 


a(x1,x2)  =  4{[c'(x1 )  -  v^x^MXj,)  -  vz(x2)]}/c02,  (1 3) 

rather  than 

o(Xi,x2)  =  ^'(x,)  ][c'(x2)]}/c02.  (14) 


In  the  next  section,  we  discuss  the  ocean  conditions  under  which  we  expect  the 
random  velocity  fluctuations  to  have  a  significant  effect  when  compared  to  the  sound- 
speed  fluctuations.  We  also  shall  briefly  survey  the  measurements  that  are  now 
available.  Finally,  we  shall  point  out  that  if  an  inverse  procedure  is  used  to  determine 
the  lateral  correlation  function,  o  ,  then  for  many  problems  it  is  not  necessary  to 
distinguish  between  the  effects  of  sound-speed  fluctuations  and  random  velocity 
fluctuations. 
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Application  to  Ocean  Problems 

Much  of  the  literature  on  ocean  scattering  of  acoustic  waves  is  devoted  to 
scattering  by  internal  waves.  In  this  case,  as  Flatte18  points  out,  the  sound-speed 
fluctuations  are  of  the  order  of  a  meter/sec  while  the  fluid  velocities  are  of  the  order  of 
centimeters/sec.  Thus  the  velocity  fluctuations  may  safely  be  neglected. 

If  one  is  interested  in  the  scattering  of  high  frequency  acoustic  radiation,  say 
20-100  kHz,  then  the  internal  waves  are  generally  no  longer  of  paramount  interest. 
Rather,  one  considers  scales  of  the  order  of  say, 10  cm  to  10  m.  In  this  case  the  sound- 
speed  fluctuations  are  of  the  order  of  perhaps  10  cm/sec  or  less,  while  the  fluid 
velocities  remain  of  the  order  of  centimeters/sec.  Hence  we  must  consider  the  velocity 
fluctuations  in  scattering  experiments,  especially  if  we  are  interested  in  the  effect  of  the 
10-cm  scale  range.  The  variability  of  the  ocean  environment  is  so  great  and  the 
measurements  so  sparse  that  it  is  dificult  at  this  stage  to  give  more  than  rough  orders 
of  magnitude  for  the  fluid  velocity  fluctuations.  We  obtained  the  above  figures  from  a 
survey  of  Refs.  1 9-24. 

Knowledge  of  the  speed  of  sound  fluctuations  in  the  10  cm  to  10  m  range  is 
greater  than  our  knowledge  of  the  velocity  fluctuations  (see,  for  example,  Refs.  25-26), 
but  still  is  hardly  satisfactory.  At  the  present  time,  we  really  do  not  have  sufficient  data 
to  make  good  predictions  of  the  intensity  and  intensity  fluctuations,  using  Eqs.(1 1)  and 
(12),  even  if  we  assume  that  sound-speed  fluctuations  are  the  dominant  effect. 
Requiring  an  independent  knowledge  of  velocity  fluctuations  is  probably  beyond  our 
present  measurement  capabilities. 

The  problem  of  finding  a  ,  when  velocity  fluctuations  are  present,  is  further 
complicated  by  the  fact  that  Eq.  (13)  shows  us  that  we  must  not  only  know  the  statistics 
of  the  velocity  field  itself,  but  must  also  know  the  function 

{c,(x1)vz(x2)}. 

That  is,  we  must  know  the  correlation  btween  the  sound-speed  fluctuations  and  the 
horizontal  range  component  of  the  velocity  field.  At  present,  we  have  virtually  no 
knowledge  of  this  function.  However,  if  the  sound-speed  fluctuations  and  the  velocity 
fluctuations  are  weakly  correlated  and  this  function  can  be  set  equal  to  zero,  then  the 
situation  is  much  better.  In  that  case  the  velocity  fluctuation  effect  is  of  order  (vz/c')2 
and  even  if  0(vz/c')  =  V3,  we  can  neglect  the  effect  of  velocity  fluctuations. 

Finally,  we  should  like  to  point  out  that  there  is  a  method  for  determining  a  which 
is  independent  of  the  relationship  between  c'  and  vz.  This  is  the  inverse  method  for 

determining  5  proposed  in  Ref.  10.  The  basic  idea  is  to  measure  the  intensity 
distribution  or  coherence  function  from  a  known  acoustic  source  and  then  to  invert  the 
solution  to  find  o.  For  example,  the  solution  for  the  coherence  function  for  a  point 
source,  in  the  absence  of  a  background  profile,  is 
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i 

{r(s,z)}  =  (I z/z2)  exp(-k2z[5(0)  -  J  5(as)da]).  (15) 

0 

Here  l2  is  a  normalization  constant  and  s  =  x1  -  x2. 

Eq.  (15)  may  be  inverted  to  find  a.  The  inversion  formula  is 

5(sx)  =  5(0)  +  d/dsx[(sx/k2z)  log({r(sx,z)}/(l2/z2))],  (16) 

with  sy  =  0.  There  is  a  similar  expression  for  sy  when  sx  =  0. 

This  method  was  used  by  Coles  and  Frehlich28  to  measure  5  for  light  propagation 
in  the  turbulent  atmosphere.  We  have  less  experience  measuring  the  coherence 
function  in  the  ocean,  but  it  appears  to  be  possible.  If  the  coherence  measurement  is 
impractical,  however,  we  may  use  instead  a  directed  Gaussian  beam  as  a  source.  In 
this  later  case,  only  the  intensity  need  be  measured  in  order  to  be  able  to  invert  the 
solution  to  find  5  .10 


Summary 

In  this  paper  we  have  shown  how  to  take  into  account  the  effect  of  a  random  fluid 
velocity  field  on  the  scattering  of  acoustic  radiation  in  the  ocean.  The  analysis  is 
subject  to  the  restrictions  (a)-(c)  given  in  the  Introduction.  It  is  then  pointed  out  that  in 
high  frequency  scattering,  where  scales  of  the  order  of  meters  become  important,  it 
may  be  necessary  to  consider  the  effect  of  velocity  fluctuations.  Finally,  the  inverse 
procedure  of  determining  o  is  discussed.  This  method  does  not  require  independent 
measurements  of  sound-speed  and  fluid  velocity  fluctuations. 
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Appendix 

We  write  Eq.  (1)  in  dimensionless  form  using  /z  as  a  characteristic  range 
coordinate  and  4-  as  a  characteristic  transverse  coordinate.  We  denote  U  as  a 
characteristic  velocity.  We  then  have 

(1/(k /z)2)  d2p/dz2  +  (1/(k/T-)2)  +  n2p 

+  2iM(1/(k/z))n2vz  dp/dz  + 

2iM(1/(k/r))n2vT»VrP  =  0,  (A1) 

where  n  =  cjc,  M  =  U/c0  and  the  independent  variables  appearing  here  are 
dimensionlized  by  ^  and  /z. 

Under  the  transformation 

p  =  u  exp(ik/zz) 

Eq.  (A1)  becomes 

(n2  -  1)u  -  2Mvzn2u  +  2i(1/k/z)  du/dz  +  (1/(k /z)2)  02u/3z2 
+  (1/(k/r)2)VT2u  +  2iM(1/k/z)n2vz  du/dz 

-  2iM(1/k/r)n2vT'7jU  =  0.  (A2) 

In  the  limit  k/z»  1 ,  the  term  (1/(k/z)2)  d2u/dz2  may  be  neglected  and  Eq.  (A2)  is 
parabolic.  However,  in  order  that  Eq.  (A2)  include  the  effects  of  index  of  refraction 
variations  and  diffraction,  we  assume  that  that  other  terms  are  of  the  same  order  as 
(1/k/z)  du/dz.  Thus,  for  example,  the  longitudinal  and  transverse  scales  are  related  by 
the  equality 

k4  =  (k*r)2  (A3) 

or  /z  =  kly2.  Furthermore,  we  demand  that  the  term  (n2  -  1  )u  is  of  the  same  order  as  the 
term  2i(1/k/z)  du/dz.  Then,  if  we  write 

(n2  -  1 )  =  2(c  -  c0)/c0  =  2n'  =  <n'2>1/2  f(x),  (A4) 

where  <n'2>1/2  is  a  measure  of  the  size  of  the  variations  in  n\  we  have 

<n'2>1/2  =  1/k/z.  (A5) 
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Eq.  (A5)  represents  a  choice  for  /z.  It  also  represents  a  restriction  on  the  size  of 
<n'2>1/2  since  we  have  already  required  that  1/k/z«  1.  This  is  the  restriction  to  weak 
scattering.  In  addition,  since  now  /j-//z=  <n'2>1/4,  we  have  the  restriction  to  forward 
scattering. 

Here,  we  also  require  that  M  «  1  and  that  M  =  m<n'2>1/2,  where  m  =  0(1).  With 
this  scaling,  the  last  two  terms  in  Eq.  (A2)  may  be  neglected.  We  find,  finally, 

2i  Bu/Bz  tV-y^u  +  2(n'  -  Mvz)  =  0,  (A6) 

which  leads  to  Eq.  (2). 

Finally,  we  note  that  the  last  term  in  Eq.  (A2)  is  the  largest  term  neglected  when 
m  =  0(1)  and  sets  the  limit  on  the  validity  of  the  approximations. 
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